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ABSTRACT: We review the properties of energetic parton propagation in hot or cold QCD 
matter, as obtained in recent works. Advances in understanding the energy loss - collisional and 
radiative - are summarized, with emphasis on the latter: it features very interesting properties 
which may help to detect the quark-gluon plasma produced in heavy ion collisions. We describe 
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two different theoretical approaches, which lead to the same radiated gluon energy spectrum. 
The case of a longitudinally expanding QCD plasma is investigated. The energy lost by a 
jet with given opening angle is calculated in view of making predictions for the suppression 
(quenching) of hard jet production. Phenomenological implications for the difference between 
hot and cold matter are discussed. Numerical estimates of the loss suggest that it may be 
significantly enhanced in hot compared to cold matter. 
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1 INTRODUCTION 

Over the past few years, a lot of work has been devoted to the propagation 
of high energy partons (jets) through hot and cold QCD matter. The jet p±- 
broadening and the gluon radiation induced by multiple scattering, together with 
the resulting radiative energy loss of the jet have been studied. These studies 
are extensions to QCD of the analogous QED problem considered long ago by 
Landau, Pomeranchuk and Migdal Jl], |j. Recent measurements ||] (reviewed in 
Q) confirm the theoretical predictions in the QED case to good accuracy. 

As in QED, coherent suppression of the radiation spectrum takes place when 
a parton propagates in a QCD medium. New and interesting predictions are 
found. When a high energy parton traverses a length L of hot or cold matter, 
the induced radiative energy loss is proportional to I? . The energy loss of a high 
energy jet in a hot QCD plasma appears to be much larger than in cold nuclear 
matter even at moderate temperatures of the plasma, T ~ 200 MeV. 

The order of magnitude of the effect in hot matter compared to the case of cold 
nuclear matter may be expected to be large enough to lead to an observable and 
remarkable signal of the production of the quark-gluon plasma (QGP). Indeed, it 
has been proposed to measure the magnitude of "jet quenching" in the transverse 
momentum spectrum of hard jets produced in heavy ion collisions, noting that 
jet quenching is the manifestation of energy loss as seen in the suppression and 
change of shape of the jet spectrum compared with hadron data. 

This review]] is organized as follows: The case of elastic parton scattering giving 

rise to the collisional energy loss, especially in hot matter, is presented in section 2. 
We concentrate on the more recent theoretical advances. References to earlier work are 

found in the quoted papers. 
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In section 3, we give the basic elements of the equations and describe the coherent 
pattern of the gluon radiative spectrum induced by multiple scattering. We 
derive the induced energy loss and the jet transverse momentum broadening 
in terms of phenomenologically significant quantities. Section 4 is devoted to 
the path integral approach, which provides another derivation of the induced 
radiative spectrum. For heavy ion collisions the case of an expanding QCD 
plasma is more realistic, and therefore in section 5 we consider the corresponding 
energy loss calculation. In section 6, we investigate the angular dependence of 
the radiative gluon spectrum. The dependence of the energy loss on the jet cone 
size is analyzed. Section 7 is devoted to estimates of the energy loss in hot QCD 
matter and in nuclear matter, and orders of magnitude are given. A noticeable 
result is that the energy loss in the case of a hot QCD medium is found to be 
quite collimated. Experimental indications are shortly reviewed. We close with 
an outlook. 

2 COLLISIONAL ENERGY LOSS IN QCD 

The electromagnetic energy loss of a charged particle in matter is a well studied 
subject |5], g. Similar mechanisms are responsible for the energy loss of a fast 
quark or gluon (jet) propagating through dense QCD matter. 

In this section we discuss the loss caused by elastic collisions of the propagating 
quark or gluon off the (light) partons forming the dense quark-gluon plasma 
(QGP). In order to understand the characteristic features we consider in some 
detail the loss of a test quark Q traversing a plasma with quarks q and gluons g 
interacting elastically as Qq — ► Qq and Qg — > Qg 0; for a review, see ||. 

The energy loss per unit length depends on the density p p of the plasma con- 
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stituents p (with momentum k) and on the differential cross section weighted 
by the energy transfer u = E — E' , where E(E') is the energy of the incoming 
(scattered) Q, 

J denotes the flux factor, q 2 the invariant (four) momentum transfer. Small 
values of q 2 dominate the collisions, 

~ C s (2 2) 

dq 2 ~ P {q 2 ) 2 ' 1 ' 

N 2 — l 

with C q = 2 C N '2 , Cg = 1 for N c colors. For a QGP in thermal and chemical 
equilibrium the densities are given by 



*N c N f m 2(N 2 - 1) 

Pi = -j2nj*~ nF< < k >> Pa = ( 27r )3 ng ( fc )' ( 2 - 3 ) 



in terms of the Fermi-Dirac (Bose-Einstein) distributions nj?(nfi). Although the 
factor cj in (^^) improves the Rutherford singularity of ( |2.2| ) , a logarithmic depen- 
dence still remains after the ^-integration, which has to be screened by medium 
effects, i.e. with a cut-off related to the Debye mass ||: — q^ ain ~ m 2 D = 0{a s T 2 ). 
Noting that Juj ~ |L (when E, » fc ~ 0(T)) one obtains @, |] 

dE 2 V- ^ M 3fc /,m 9max ^OL 2 T 2 ( ^ Nf\ . cE , nA . 

-^ = ^T.C P J — Pp{ k) In S- * (l + J ) In — , (2.4) 

with c a numerical constant of 0(1), and N c = 3. The strong coupling constant 
may be evaluated at the scale a s (T) for high temperature T. 

Because of the T 2 dependence of (2.4) it has been pointed out by Bjorken 
|7j that the collisional loss is proportional to y/e, i.e. the square root of the 
QGP energy density, which in leading order in the coupling constant is given by 
e = 8vr 2 T 4 (l + 21JV//32)/15 0. 
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A proper and consistent treatment of the screening effects of the plasma in 
the low (soft) exchange momentum region of the collisions is indeed possible in 
the thermal field-theoretic framework Q, using resummed perturbation theory 
at high temperature. This method has been developed by Braaten and Pisarski 
[jll], 12 1 and it allows one to calculate the hard thermal loop (HTL) corrections 
to the propagator of the exchanged gluon in the Qq — ► Qq, and the Qg — > 
Qg processes. The quantum field-theoretic calculation of the energy loss of a 
quark requires the evaluation of the discontinuity of the self-energy diagrams e.g. 
illustrated in Figure |]. For the soft momentum exchange (with momentum less 
than q c = 0(g l/2 T)) the HTL gluon propagator contributes, whereas for the 
hard momentum exchange (q ~ T) the tree-level elastic scattering (Figure [Tjb) 
contributes [13]. 

The momentum cut-off q c drops out in the sum of soft and hard contributions. 
It is important to note that Landau-damping effects, because of the negative q 2 
values in Q2.1| ), screen successfully the low q 2 region leading to a well defined 
result for dE/dz (at least to leading order in the coupling constant). 

As an example the result is illustrated in Figure 0, where the energy loss of a 
charm quark is shown by the dashed curve using parameters characteristic for a 
thermalized QGP as expected in ultrarelativistic heavy ion collisions. 

For light quarks the collisional energy loss for a jet propagating in a hot medium 



of T = 0.25 GeV amounts to 0.2 - 0.3 GeV/fm g, gj, in agreement with the 
estimates of |?| as shown in Figure (dotted curve). For a gluon jet the loss is 
predicted to be larger by the color factor 2N 2 /{N 2 — 1) = 9/4 than for the quark 
jet. 

Since the QGP expected at RHIC and LHC is likely to be out of chemical 
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equilibrium it is necessary to investigate the energy loss in this case [Tq] . 
Indeed, even away from chemical equilibrium, dynamical screening remains op- 
erational within the HTL-resummed perturbation theory. More explicitly, the 
collisional energy loss for a heavy quark (mass M) propagating through a QGP 
parametrized in terms of the distribution functions X q np and X g riB, respectively, 
where X q ^ g are the fugacity factors describing chemical non-equilibrium, becomes 



dE 

— — = 2a^m n In 
dz 9 



0.020 V ^ 2 x " N f^ 12 ^+ 2 \ N f) 



(2.5) 



This expression [16] is valid for energetic quarks with E 3> M /T and contains 



for X q = X g = 1 the original result of [13]. The screening mass parameter is 



™2 



in 
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ATra s (X g + X q N f /6)T 2 /3. (2.6) 



For comparison the solid curve in Figure |2| shows the loss for the interesting 
case of the "early plasma phase" which is dominated by gluons (X g = 1, X q = 0), 
where the loss is exclusively due to elastic Qg — > Qg scattering mediated by gluon 
exchange. 

In summary, even when the partons propagating in hot matter have a large 
momentum, the collisional energy loss per unit length turns out to be less than 
0(1) GeV/fm when reasonable values for a s and T are taken. This estimate 
may be compared with the value for the hadronic string tension, k — 1 GeV/fm, 
which measures the slowing down of a high momentum quark in (cold) nuclear 
matter [01. 



3 RADIATIVE ENERGY LOSS IN QCD 



3.1 Model, basic parameters and equations 

We imagine a very energetic quark of energy E propagating through a QCD 
medium of finite length L. Multiple scattering of this projectile in the medium 
induces gluon radiation, which gives rise to the quark energy loss. 

The main assumption [18, 19] is that the scattering centers are static and 



uncorrelated (in the spirit of the Glauber picture). We thus focus on purely 
radiative processes since the collisional energy loss vanishes in the case of static 
centers. 

We define a normalized quark- "particle" cross-section 

V{Q 2 ) = (3-7) 

7T<T aQ z 

where Q is the 2-dimensional transverse momentum transfer scaled by an appro- 
priate scale : 

<?-*, 

and a = J dQ . (3.8) 

In the case of a hot QCD plasma, the "particle" is a quark or gluon and it is a 
nucleon in the case of cold matter. da/d 2 Q depends only on Q, as it is usually 
assumed for diffractive kinematics with very large incident energy. The scale u 
characteristic of the medium is conveniently taken as the Debye screening mass 
in the hot case and as a typical momentum transfer in a quark-nucleon collision. 
The condition that the independent scattering picture be valid may be expressed 
as : 

a' 1 < A, (3.9) 

where A is the parton mean free path in the medium A = 1 /pa, and p is the 
density of the medium. We assume that a large number of scatterings takes 



place, that is 



L> A. 



(3.10) 



Successive scatterings being independent, the parton propagation is "time- 
ordered" and time-ordered perturbation theory is the natural framework to cal- 
culate the radiation amplitude. Let us give a sketch of the basic equations, 



N N 



(3.11) 



referring the reader to [20, 21] for further details. We may number the scattering 
centers depending on the interaction time and write for the radiation spectrum 
induced by N scatterings 
dl a v 

2 / dk± [ 

\ i=l 3=1 

where is an effective current for the gluon emission induced by center i. It 
includes color factors consistent with the overall normalization to the elastic scat- 
tering cross section. The phase ifi is associated to time U (longitudinal coordinate 
zi) by ifi = tikj_/iv. 

The brackets indicate averaging - over momentum transfers and over Z{ - for 

which a simplified model is 

r N-l i a /A\ r N 

(( - )) <*/ n ^ e -p (-t) 7n^ v ^ ■■■ ) • ( 3 - 12 ) 



where = — 2^. We rewrite (3.11) as 
dl 



UJ- 



duo 



2tt 2 



dkj_ 



I N N 

2 ^E E 

\ i=lj"=i+l 



1 + 



N 



E J ~ 

i=l 



(3.13) 



which allows one to exhibit the so-called factorization contribution: the second 
term in ( 3,13j ) which corresponds to the limit of vanishing phases. This contri- 



bution is equivalent to the radiation spectrum induced by a single scattering of 

N 

momentum transfer q± tot = Qi- It nas & t most a weak logarithmic medium 

i=i 
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dependence We concentrate in the following on the medium- induced radi- 
ation spectrum and drop the factorization term. In the limit of large--/V c and 
replacing sums over i and j by integrals, i.e. taking the sum over scatterings to 
be arbitrary in number, the following expression for the spectrum per unit length 
of the medium is obtained, 



dl a s N c f , rL—A dzi 



dujdz 2tt 2 L Jo Jq A 



dA dU 



2Re ]T J x ■ J n+2 



n=0 



ex P y K l^ u e—j ~ l 



n+l 



8 A-^A W ), (3.14) 



m=l 



with U = k±/ fi and k = A/U 2 /2uj. In the soft gluon limit the rescaled emission 
current is given by 



^=(^-77^d- (3-15) 



Ui Uj + Q 
Pi (Ui + Qi) 2 

Eq. ( |3.14| ) exhibits the interference nature and coherent pattern of the spectrum. 



The phase factor as it appears here may be understood in terms of formation 
time arguments which will be discussed heuristically in the next section. It can 
be shown moreover [21] that (3.14) leads to a simple structure of the spectrum 



per unit length which will be the starting point of section 3.4: 

dl a s N c [L dA rL-A dzi r ? - . * =0 

K 

where fBom(U) is the Born amplitude defined as 



(3.16) 



fBorn(U) = J dQ x V(Q\) J x , (3.17) 
and f(U, A) is the evolved amplitude which satisfies a Bethe-Salpeter type equa- 



tion. Subtracting in ( p. 16 ) the contribution for k = corresponds in ( p. 14 ) 



to subtracting the zero phase contribution. The generic structure of (3.16) is 
illustrated in Figure ||. The reader may question the one-gluon exchange approx- 
imation as shown in this Figure |3[ In fact, the scale of the coupling for each 
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individual scattering is set by the accumulated overall transverse momentum. 
This justifies the perturbative treatment. 



3.2 Heuristic discussion 

In the following discussion we neglect logarithmic and numerical factors of 0(1), 
but keep all relevant parameters. 

The semi-quantitative argument which allows one to understand the coherent 
pattern of the induced gluon radiation is the following. One defines the formation 
time of the radiation, 



£form — ,2 ' (3.18) 
K I 



where u and k± are the gluon energy and transverse momentum (with to 3> k± 
and the typical k± ~ fj,). When t{ OTm S> A, radiation takes place in a coherent 
fashion with many scattering centers acting as a single one. Let us introduce 
the coherence time (length) Z co h which plays an important role in the following 
considerations. It is associated with the formation time of a gluon radiated by a 
group of scattering centers which acts as one source of radiation, 

^coh — 77 2 \ ' (3.19) 



with 



(kl)i coh = l -f f = iW, (3.20) 
assuming a random walk expression for the accumulated gluon transverse mo- 
mentum. One derives the estimate 



/ coh -W4^, (3.21) 



so that the number of coherent scatterings becomes 
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where the energy parameter -E"lpm = -fy/ 2 is introduced , in analogy with the 
QED Landau-Pomeranchuk-Migdal (LPM) phenomenon. 

For small u < .Elpm; incoherent radiation takes place on L/A scattering cen- 



ters. Using the soft oj limit for the single scattering spectrum [23] 



udl a s 
— — ~ — N c , 

dw 7T 



(3.23) 



the differential energy spectrum per unit length in the so-called Bethe-Heitler 
(BH) regime for incoherent radiation is derived, 



ujdl 



dw dz 



1 udl 



BH 



L dw 



- vr C A' 



(3.24) 



with / co h < A and u < o>bh = -E'lpm- 

The interesting regime of coherent radiation ( LPM regime ) is defined by A < 
/ coh < L (N coh > 1), i.e. 



Elpm = w B h < w < min {uj {£LCt ,E}, 



(3.25) 



with Wf a ct ~ \L 2 . Since the A?coh groups are acting as effective single scattering 



centers the energy spectrum is estimated as 



ujdl 



duo dz 



1 udl 



LPM 



/ co h dw 



^coh 



a s , r 1 a, , r / u? 1 

7T f co h 7T \l \ UJ 



(3.26) 



Comparing (pT26| ) with ( gjg ) we find a suppression factor given by ^/E'lpm/w. 
For / C oh > ^) i- e - when 



^ > ^fact = -^LPM 



(3.27) 



effectively only one scattering is active ( factorization regime ) . and correspond- 
ingly for w fa ct < u < E, 



uodl 



dw dz 



fact 



— N c \. (3.28) 

7T L 
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The expressions for the spectrum per unit length in the different regimes 
(Eqs.(|3~24l), (j3T26j) and ( ^28| )) hold for a medium of finite length L < L cr , where 



L cr = Ay -E/-Elpm, (3.29) 

as derived from the condition that u;f act < E (correspondingly, the condition 
E > E cr = Elpm^/A) 2 has to be satisfied). A discussion of the radiation 



spectrum can be found in [21]. 



In order to obtain the energy loss per unit distance —dE/dz one integrates the 
gluon spectrum over u, with < to < E. In addition to a medium independent 
contribution proportional to —N c j^ (the factorization contribution), we obtain 



from (3.26|) the medium induced (LPM) loss, proportional to the size of the 



medium and given by 



12 2 

- dE/dz ~ —N C J ^-u iact ~ ^N c ^-L, (3.30) 

7T y A 7T A 

for L < L cr . Integrating over z leads to the total loss growing as I? . For L > L cr 
(i.e. E < E cr ), the size does not affect the loss per unit length, 



dE/dz ~ *N c Me = ^ ^^E^E, (3.31) 

7T V A 7T A 



i.e. a dependence proportional to yE is obtained, which is familiar from the 
QED-coherent LPM suppression Q. 
In Figure || the energy loss 

pL 

-AE= --z-dz, (3.32) 
Jo dz 

for the induced and the factorization cases, is shown as a function of L (iV c is 
taken to be 1). 

Using the random walk expression for the accumulated transverse momentum 
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of the gluon due to successive scatterings in the medium of size L 



(k 2 ± ) L ~ fi 2 L/ A, 



(3.33) 



and inserting this relation in ( 3.30| ), 



one obtains 



dE/dz ~ — N c (k 2 ± ) L , 



(3.34) 



a relation between the induced energy loss and the jet broadening, which is 
independent of the details of the interaction. 

3.3 Jet p ^-broadening 

On the way to deriving the gluon radiative spectrum, let us start with the clas- 
sical diffusion equation satisfied by the transverse momentum distribution of a 
high energy parton which encounters multiple scattering in a medium. Suppose 
the parton is produced in a hard collision with an initial transverse momentum 
distribution fo{U 2 ) ; U is the dimensionless transverse momentum U = p//i and 



Neglecting the transverse momentum given to the parton by induced gluon 
emission, one can derive a kinetic equation for the transverse momentum distri- 



In terms of the variable t = z/\r with A# the mean free path for a parton of 
color representation R, one finds the following gain-loss equation 



JdUf (U 2 ) = l. 



bution f(U 2 ,z) after a distance z in the medium [p4|]. 




= + 




dt 




(3.35) 



with 




(3.36) 
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Defining f(B 2 ,t) as 

f(B 2 ,t) = J dU e~ iM f(U 2 ,t), (3.37) 



and 



we find 



with 



V(B 2 ) = J dQ e~ i3 $ V(Q 2 ), (3.38) 



df(B 2 ,t) 1 d2 ^ d2 w, d2 



dt 4 



-B 2 v(B 2 )f(B 2 ,t), (3.39) 



v(B 2 ) = ±(l-V(B 2 )). (3.40) 



It is possible to define a characteristic width of the distributions f(U ,t) which 



is found to be [g4| : 



Plw = v(\ R /L). (3.41) 

The linear growth with L is expected and is used to discuss pj_-broadening of high 
energy partons in nuclei. The coefficient q = plays the role of a transport 
coefficient as encountered in diffusion equations. (3.41) is valid for hot and cold 
QCD media. 

3-4 The radiative gluon energy spectrum and induced energy loss 

Let us now turn to the gluon spectrum 0. We shall concentrate on a quark 



jet. The general case is given in [25]. As sketched in section 3.1 and derived 



in |2(|, 21 , ^] , the spectrum for the radiated gluon is calculated in terms of the 
interference term between the quark-gluon amplitude at time t and the complex 
conjugate Born amplitude. For simplicity, we restrict ourselves here to the case 
where the quark enters the medium from outside. (An additional term is needed 



"For related discussions of gluon bremsstrahlung in dense matter see also [E6l G 
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in the case when the quark is produced via a hard scattering at t = in the 
medium). We denote by f(U,V,t) the quark gluon amplitude at time t. U is 
the scaled gluon momentum U = £ and V — U the scaled quark momentum as 
illustrated in Figure ||. 

To account for the gluon polarization, / is a 2-dimensional vector which will be 
implied hereafter. The dependence on U and V is actually only in the combination 
U — xV with x = k/p. The amplitude f(U,V,t) satisfies the initial condition 
f(U,V,0) = f Bo rn(U,V), where f B orn is the Born amplitude to be described 
shortly. 

The induced gluon spectrum is written as : 

-2Re I dU < I dt 2 I dh 



du) dz ir 2 L 

P<y^f(u-xV,t 2 -t 1 



o Jo 



(3.42) 



The various terms in ( 3.42; ) have simple interpretations. The as J> F is the coupling 
of a gluon to a quark. The 1/L comes because we calculate the spectrum per 
unit length of the medium. The factor 2TT^f(U — X V, t 2 — t\)padt\ is the number 
of scatterers in the medium, padti, times the amplitude with gluon emission 
at ti, evolved in time up to t 2 , the time of emission in the complex conjugate 
amplitude. The factor fB orn (U — xV)padt 2 gives the number of scatterers 
times gluon emission in the complex conjugate Born amplitude. The subtraction 
of the value of the integrals at u = oo eliminates the medium independent zero- 
phase contribution. Eq. ( 3.42] ) may be simplified using t = (^FA^) r. Defining 
ro = 2 %^j, we obtain 

ujdl a s N c 



dw dz 7r 2 A 



-Re J dQ 



dr (l - Q f(U - xV,t) ■ f* Born (U - xV)Y " • (3.43) 
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Due to the specific dependence of / and fsom on U and V, it is possible to 
express them in terms of a single impact parameter as : 

fBorn(U -xV)= J -^L e iS <°~^ f Born (B) , (3.44) 

allowing us to obtain the following expression for the spectrum in impact param- 
eter space : 

£5-^*/i{jr*( i -s)^-«-*r- ™ 

The generic diagram appears in Figure ||. The complete list of diagrams describ- 
ing the Born amplitude is shown in Figure ||. Graphs a-c correspond to inelastic 
reactions with the medium while graphs d-g correspond to forward scattering 
in the medium. For terms a-c there are corresponding inelastic reactions in the 
complex conjugate amplitude. In the approximation that the forward elastic am- 
plitude for quark scattering off particles in the medium is purely imaginary, the 
elastic and inelastic terms are proportional to V(Q 2 ). The color factors and the 



expression of each graph contribution are derived in [25]. 



The quark-gluon amplitude f(U, V, t) obeys an integral evolution equation 



derived in [2C, 21, 25 1. In impact parameter space and in the small- x limit, this 



equation takes the simple form 

^-f(B, t) = in V| f(B, t) - 2(1 - V(B))f(B, r) (3.46) 

with k = (^uTj an< ^ f(B,0) = fBom{B). This equation is a Schrbdinger- 
type evolution equation for the propagation of the quark-gluon system in a QCD 
medium. Comparing ( p. 46 ) to ( p.39| ) is instructive. The term proportional to k in 



(3.46) is clearly of quantum origin and is associated to the phase of the amplitude 
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whereas ( |3.39| ) is a classical diffusion equation. The contributions entering the 
expression for the spectrum (3.43) are depicted in Figure [/]. 



So long as v(B 2 ) = 4(1 — V(B)/B 2 can be treated as a constant, solving 
( |3,46| ) proceeds in analogy with that of the 2-dimensional harmonic oscillator 
with imaginary frequency. We expect that the behavior of v(B 2 ) is close in 
general to the Coulomb potential case i.e. ~ £n(l/B 2 ) at small B 2 . The solution 



of ( 3.46 ) to logarithmic accuracy is worked out in [21, f25| . 

In the case where the (quark) jet is produced in matter, one finds for the gluon 
spectrum 

udl 2a s C F x 2 

[1 - x + — \ ln|cos(w To)| , (3.47) 



dwdz irL 2 
from which the energy loss per unit length is derived, 

dE r 00 uj dl 



dz Jo doo dz' 



(3.48) 



i.e. 

,2 



dE a s N c n*L _ x 

Tz = — — u(T ° } - (3 ' 49) 



Notice the remarkable relation (cf. Eq. ( 3.34 )) 



dE u s N c 2 
- Iz- = (3-50) 

between energy loss and jet p± -broadening |24j| . 
4 PATH INTEGRAL APPROACH 

This section is devoted to presenting in some detail a different approach to the 
LPM effect in QCD, based on the so-called light-cone path integral technique 



for multiple scattering developed in [28]. We refer the reader for derivations and 



further details to [28] - [p3j] . Let us give here the essential features of the method 
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by discussing in scalar QED the formalism for an induced transition a — * be. The 
interaction Lagrangian is Li n t = A ipliplipa + h.c . It is assumed that the decay 
a — ► be does not take place in the vacuum. We shall indicate later the proper 
treatment for realistic QED and QCD. 

4-1 Derivation of the basic formulas 

The /S-matrix element for the a — > be transition in an external potential reads 

(bc\S\a) = i [ dtdfX(z) ifj* b (t,r)4>*(t,r)Mt,7), (4.51) 



where ipi are the wave-functions, and \{z) adiabatically vanishes at \z\ — > 00 . 
Let us consider the case of a static external potential. Then we can write tpi as 

^i(t,r) = exp[-iEi{t - z)]</>i(z,p), (4.52) 

where p is the transverse coordinate, and the function 0j describes the evolution 
of the ipi on the light-cone t — z = const. At high energies Ei 3> rrii, after 
substituting ( 4.52| ) into the Klein-Gordon equation, one can obtain for <fii the 



two-dimensional Schrodinger equation 



ij£ = Hi(z)^, (4.53) 



H ^=-^XU + ^ ulp - z)+n i,' (4 - 54> 

where fii = E{, e^ is the electric charge, and U is the potential of the target. 
Consequently, the values of the fa at the p-planes z = Z2 and z = Z\ are related 

by 

<f>i(z2,p2) = J dp 1 Ki(p 2 ,z 2 \pi,zi)(j)i(zi,p 1 ), (4.55) 
where Ki(p2, Zi\p\, z\) is the Green function of the Hamiltonian (4.54). Let us 



introduce the two p-planes located at large distances in front of (z = Zj) and 
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behind (z = Zf) the target. Then, using the convolution relation ( 4.55| ) one can 



express the incoming and outgoing wave-functions in terms of their asymptotic 
plane- waves at Z{ and Zf, respectively. As we shall see below, this representation 
turns out to be very convenient for the evaluation of the LPM effect. It is one of 
the key points of the light-cone path integral approach. 
The differential cross section can be written as 

d 5 a 



dxd%dq c (27r)' 



-Re J dpidp 2 J dz 1 dz 2 gF(z 1 ,p 1 )F*(z 2 ,P2), (4.56) 

Zl<Z 2 

where F(z,p) = 4>l(z, p)fa(z, p)4> a {z, p), % jC are the transverse momenta, x = 
Eb/E a , and g = X(zi)\(z2)/[1Qttx(1 — x)E%] is the vertex factor. Expressing fa in 
terms of the asymptotic plane- waves, ( |4.56| ) may be represented diagrammatically 
by the graph of Figure [8|a. We depict Ki{K*) by — > {*—). The dotted lines show 
the transverse density matrix for the initial particle a at z = Zi, and the complex 
conjugate transverse density matrices for the final particles b, c at z = Zf. For 
the spectra integrated over transverse momenta, the relation 

J dp 2 K(p 2 , z 2 \pi,z 1 )K*(p2, z 2 \p[,z 1 ) = 5(pi - p[) (4.57) 

allows one to transform the graph of Figure ||a into the ones of Figure ||b and 
Figure ||c for the q c - and <fb jC -integrated spectra, respectively. 

Let us discuss now the a — ► be transition for a random potential of an amor- 
phous target using the representations of Figure ||. In this case one should perform 
averaging of the transition cross section over the states of the target. We cannot 
evaluate analytically the diagrams of Figure |8] for a given state of the target. The 
basic idea of the approach of [29, 30, 34] is to represent all the propagators in the 
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path integral form 

Pi(dp/dz) 2 



Ki(p 2 , z 2 \pizi) = J Dpexpli J dz 



eiU{p,z) 



im 2 (z 2 - Z\, 



2pi J ' 
(4.58) 



and to perform averaging over the target states before the integration over the 
trajectories. It is then remarkable that for the diagrams of Figure ^o,c (and for 
Figure |8)a if b or c has zero charge, say, for the e — > je' transition) a considerable 
part of work on the path integration can be done analytically. 

Below we consider the ^-integrated spectrum. Let z\ and z 2 be the longi- 
tudinal coordinates of the left and right vertices of the diagram of Figure |8|b. 
Taking advantage of the convolution relations (we omit the transverse variables) 
K h (zf\z\) = K h (zf\z 2 ) <8> K b (z2\zi) , K*(z 2 \zi) = K*(z 2 \zi) ®K*(zi\zi), we can 
divide the diagram of Figure |8|b into the initial and final state interaction two- 
body parts (we denote them by S c and Sy and the three-body part (we denote 
it by M) located between them. The factors Si and M are given in terms of the 
Green functions ( [4.58 ). 



Let us first consider the factor Sf after averaging over the states of the target, 
the phase factor 

exp{iei J dz[U(p,z) -U(p',z)}} (4.59) 

can be viewed as the Glauber factor for the ii pair. Neglecting the correlations in 
the positions of the medium constituents one can obtain for the averaged value 
of this phase factor (we denote it by <J>j({/?}, {/?'})) 



iP'}) = ex P 



dz(Jii(\p{z) -p'{z)\)n{z) 



(4.60) 



where {p} and {p'} are the trajectories for K{ and K*, respectively, ag is the 
dipole cross section of the interaction with the medium constituent of the ii pair, 



22 

and n{z) is the number density of the target. Then the S{ is given by the path 
integral / DpDp' exp[iSi({p},{p'})] with the action 

Si({p},{p'}) = \J dz{p t [(dp/dzf - (dp'/dzf]+iad\p(z) - p'(z)\)n(z)}. 

(4.61) 

It is important that the interaction term in ( |4 . 6 1|) depends only on the relative 
distance between trajectories. This fact allows one to carry out analytically the 



path integration and obtain a simple formula [28] 

Si(p2,P2, z 2 \pi,p{,zi) = Ki(p 2 , z 2 \pi, Zi)K!-*(p2, z 2 \p[,z 1 )$ i ({p s }, {p' s }), (4.62) 
where 

m(P2 ~ Pi) 2 imf(z 2 -zi) 



Kf(p 2 ,z 2 \pi,z 1 ) = — exp 

2m(z 2 - z\) 



(4.63) 



2(^2 - z 1 ) 2pi 
is the Green function for U = 0, {p s } and {p' s } denote the straight line trajectories 



between pi j2 and p{ 2 , respectively. Expression (4.62) can be obtained by dividing 
the z-interval into steps of small width, and taking the multiple integral step by 
step [§§. 

The factor M can be treated similarly. The corresponding Glauber factor 
contains the three-body cross section Oabc depending on the relative transverse 
separations = Pb — p c and t Sc = p a — p c (here and below we view the particle 
a for a complex conjugate propagator as antiparticle a). The path integrals may 
be performed analytically leading to 

M(p 2 ,P2,z 2 \pi,p[,zi) = K°(R 2 ,z 2 \Ri,Zi)K°*(p2,z 2 \pi,z 1 )K bc (p2 - fr, z 2 \0, zi), 

(4.64) 

where R\ = p\ , R 2 = xp 2 + (1 — x)p' 2 are the initial and final coordinates of the 
centre-of-mass of the be pair, respectively. The Green function Kb c is given by a 
path integral on t% c and describes the evolution of the bed system. 
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Having fl4~6p and ( jOjgj ), one can represent the spectrum for the a — » 6c 
transition for the diagram of Figure ||b for q a = in the form [34] 

d 3 I 2 



dxdq b 



(2tt)2 



Re y dr;, exp(— iq b ■ f b ) J dz\ 



J dz 2 g <&b(rb, z 2 )K bc (T b , z 2 \0, zi)<£> a {r ai Zt), (4.65) 
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where 



$a(T a ,zi) = exp 



$b(n,z 2 ) = exp 



Paa(T a ) 



«1 



dzn(z) 



dzn{z) 



(4.66) 



are the values of the absorption factors for the parallel trajectories, and r a = xTb- 



The potential for the Green function Kb c entering (4.64) and (4.65) should be 
evaluated for parallel trajectories of the centre-of-mass of the be pair and a. The 
resulting Hamiltonian for Kb c is given by 



H, 



be 



1 / d 



in(z)aa bc (u c ,Tg C ) J_ 
2 + V 



(4.67) 



2fi bc \df bc , 

where \i bc = E a x(l — x) is the reduced Schrodinger mass, fa C = xfb c — T a ; the 
formation length is Lj = 2E a x(l — x)/[m b (l — x) + m^x — m^x(l — x)\. 



For numerical calculations it is convenient to represent ( 4.65 ) in another form 
in which (for a target occupying the region < z < L) the z-integration is 
dominated by the region \z\ < max(Lj,L). Let us rewrite the integrand of the 
^-integral in ( 4.65| ) as 



9 {$b(n,Z 2 )[K bc (T b ,Z 2 \0,Z 1 ) - K$ c (T b ,Z2\0,Zi)]$ a (T a ,Zl) 

+ [Mn,Z 2 ) - l]i^ c (f b ,Z 2 |0,Zi)[$ a (f a ,Zi) - 1] +K$ c (T b ,Z 2 \0,Zl)[<f> a (f a ,Z 1 ) - 1] 

+ [Mn, z 2 ) - l]Kl(T b , z 2 \0, z x ) + Kl(n, z 2 \0, Zl )}. (4.68) 
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It is evident that in (4.65) the contribution associated with the first two terms of 



( f4,68| ) will be dominated by the region \z% t2 \ < max(Lj, L), and can be evaluated 



neglecting the z-dependence of X(z). However, for the last three terms in (4,68) 



it is not the case. Taking z% = — oo, Zf = do, and using X(z) which exponentially 
vanishes at \z\ — > oo one can show that the term oc K bc does not contribute to 
the spectrum. It is not surprising since this term corresponds to the transition 
in vacuum. The contribution of the other two terms can be written in terms of 
the be Fock component of the light-cone wave-function of the particle a, f] 
The final expression for the spectrum is given by 

Zf Zr 

Si 2 



dxdqh (2ir) 



—Re J dr b exp(-iq b • r b ) J dz x J dz 2 

Zi 21 

• g{®b(n,z 2 )[K bc (T b ,z 2 \0,z 1 ) - K bc (T b ,Z 2 \0,Z 1 )]$ a (Ta,Z 1 ) 

+ [Mn,z 2 ) - i]K b \(T b ,z 2 %z 1 )[<s> a {T a ,z 1 ) - i]} 

- 7^A2 / dfd ^' e M~iQb ■ T)V b a c *(x, t - r)^ b a c (x, r) 

■ [$ b (T,Z i ) + $a{xT t ,Zf)-2], (4.69) 

where one can take z% = — oo, zt = oo. 

Integrating over q b one obtains from ( 4.69| ) the x-spectrum 



z f z f 

dJ_ 

dx 



2Re / dzi / dz 2 g K bc (p 2 , z 2 \pi, z{) - K$ c (p 2 , z 2 \pi, z{) _ _ , (4.70) 

J J L J Pl=P2=0 
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which was derived in [29] using the unitarity connection between the probability 
of the a — > be transition and the radiative correction to the a — ► a transition. 

The spectrum Q4.70|) can be represented in another form demonstrating a close 
connection between the LPM suppression and Glauber absorption. Treating the 

3 These terms, which vanish after integrating over qt, have been missed in 
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second term of the Hamiltonian (4.67) as a perturbation, one obtains 

Z'2 

K bc (0,z 2 \0, Zl ) = Kl(0,z 2 \0,zi)+ JdtJdpKl(0,z 2 \p,(;)v(ti,p) 

Zl 

Z2 Z2 

■K^(p^\0, Zl ) + J ' d& J d6 J ' dp l dp 2 Kl{^z 2 \p2^2)v{^P2) 



Kbc(P2,&\pi,£l)v{£l,px)K% c {px,Z,i\ti,Zi), 



(4.71) 



where v(z, p) = —n(z)<Jabc(p, xp)/2. Taking advantage of ( 4.71 ) one can represent 
(4.70) in the form 

dl _ dI BH + dl abs 
dx dx dx ' 

= T J dp\* b a c (x,p)\ 2 a- abc (p,xp), (4.73) 

L L 

d jabs 



dx 



--Re J dzin(zi) J dz 2 n(z 2 ) J dp^ b a c *(x, p)a abc (p,xp) ^(x, p, z l7 z 2 ), 

Zl 

(4.74) 



where T = Jq dzn(z) |30|], and 



$(x,p,zi,z 2 ) = I dp'K bc (p,z 2 \p\zi)wl c (x,p')aabc(p',xp') 



(4.75) 



is the solution of the Schrddinger equation with the boundary condition 



<S>(x,p,Zl,Zl) = y b a C (x,p)a a bc(P,Xp) 



(4.76) 



The first term in ( 4.72| ) corresponds to the impulse approximation. It dominates 
the cross section in the low-density limit (the Bethe-Heitler regime). The second 
term describes absorption effects responsible for the LPM suppression. 

For a sufficiently thin target the absorptive correction ( [4.74 ) can be evaluated 
neglecting the transverse motion in the dbc system inside the target (it corre- 
sponds to neglecting the kinetic term in ( [4.67| )). Then, the Green function takes 
a simple eikonal form 

Vabc(pi,xpi) f Z2 



K bc (p2, z 2 \pi,z 1 ) 5(p 2 - pi) exp 



dzn(z) 



(4.77) 
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and (gjg)-(gg) give 

dl fr 



dp\^ b a c (x,p)\ 2 T^ c (p^P), (4-78) 



dx 

with r^(f bc , f ac ) = {1 - exp[-Ta a bc(n c ,T ac )/2]}. This is the "frozen-size" ap- 
proximation corresponding to the factorization regime discussed in section 3. 

The above analysis is performed for the particle a incident on a target from 
outside. If it is produced in a hard reaction inside a medium one should replace in 
equations ( [4.65|) , ( |4.69 ) and ( 4.7C| ) z- L by the coordinate of the production point, 



and in ( [4.69D the factor [$ a — 1] should be replaced by <£ a . Note that due to the 
infinite time required for the formation of a light-cone wave- function ^ b a c , ( ^72|) 
does not hold in this case. 

4-2 Generalization to the realistic QED Lagrangian 

The generalization of the above analysis to the realistic QED and QCD La- 
grangian is simple. Let us consider first the e — > e'j transition in QED. The 



5-matrix element can be obtained by replacing in (4.51) A by eu e >^ u e u u e where 
e u is the photon polarization vector, and u e i,u e are the electron spinors in which 
the transverse momenta should be regarded as operators acting on the corre- 
sponding wave- functions. Since the photon does not interact with the target, one 
has crg 7e ' (r 7e ' , r ge ' ) = cr ( | fg e ' | ) , where a is the dipole cross section for the e + e~ 
pair. In terms of the electron-atom differential cross section it reads 



°iP) = IJ d <l\ l ~ e MiQ- P)}^2- (4-79) 



'dq 2 

The dipole cross section vanishes as p — > 0, and one can write it as cr(p) = C(p)p 2 , 



where C(p) has a smooth logarithmic dependence at small p [28, 36]. 



In the Bethe-Heitler regime the radiation rate is dominated by Tg e ' < l/m e ; for 
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the case of strong LPM suppression the typical values of rg e ' are even smaller. One 
can approximate the Hamiltonian ( [4.67 ) by the harmonic oscillator Hamiltonian, 
and obtain from ( |4.70| ) the radiation rate per unit length. In an infinite medium, 
in the regime of strong LPM suppression, the radiation rate per unit length takes 
the form 



dl _ a[4 -4x + 2x 2 } I C( Peff x) 



dxdz 2ir y 2x(l — x)E e 

The value of p e ff can be estimated as p e ff ~ (2Lj//i 7e /) 1 / 2 where the formation 
length L'j is the typical value of \z% — Z\\ in ( 4.7C| ). One can see that for strong 
suppression p e ff becomes much smaller than l/m e x, the characteristic transverse 
size in the Bethe-Heitler regime. For this reason the spectrum for strong sup- 
pression ( |4.80| ) is insensitive to the electron mass [35]. Note that the oscillator 
approximation is equivalent to the Fokker-Planck approximation in momentum 
representation used in Migdal's analysis Q. This fact is not surprising since 
within logarithmic accuracy <r(p) oc p 2 leads to a Gaussian diffusion of the elec- 



tron in transverse momentum space [28|. This feature underlies the relationship 
given in section 3 relating the energy loss and pj_-broadening in QCD. 

For an accurate numerical evaluation of the LPM effect it is convenient to use 



the form given by ( 4.72 )-( 4.74 ). In |3C , 33] it is used for the analysis of the recent 
data on bremsstrahlung from high energy electrons taken by the E-146 SLAC 
collaboration Q. Excellent agreement (at the level of the radiative corrections) 
with the data is found. 



4-3 Induced gluon emission in QCD 

Let us now discuss the induced gluon emission from a fast quark in QCD. At the 
level of the radiation cross section, involving the sum over states of the medium, 
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one can formulate the theory similarly to the case of QED. The path integral 
representations for the diagrams of Figure |8] can be written by introducing into 
the vacuum path integral formulas the Glauber factors for propagation of the 
color neutral partonic systems (consisting of the partons from the amplitude 
and complex conjugate one). The quark trajectory for the complex conjugate 
amplitude can be regarded as that of an antiquark with negative kinetic and mass 
terms. It follows from the relation — T* = T q (here T Qtq are the color generators 
for a quark and an antiquark). The qq, gg, qqg configurations which can appear 
in the graphs like those of Figure ||b,c, have only one color singlet state, and the 
diffraction operator has only diagonal matrix elements involving the two-body 
cross sections a qq (p),a gg (p) = \a qq -(p), and the three-body one cr gqq (p gq , p qq ) = 
\[°qq{\Pgq\) + a qq(\Pgq\)]-\°qq(\Pqq\) © • Thus, the spectra integrated over quark 
or /and gluon transverse momenta can be evaluated similarly to the above case 
of the a — ► be transition in QED, with all the particles now being charged. 

For the x-spectrum ( |4.70| ) p qq = xp gq , and the three-body cross section takes 
the form a gqq (p gq , xp gq ) = \[a qq (p) + a qq ((l - x)p)\ - \a qq (xp), where p = 
\p gq \- Similarly to QED, one can estimate the spectrum using the oscillator 
parametrization a gqq (p,xp) Cs(x)p 2 , where Cz(x) = |{9[1 + (1 — x) 2 \ — 
x 2 }C 2 (p e ff), C 2 (p e ff) = a qq {peff)/plff Here p eff is the typical size of the 
qqg system dominating the radiation rate, which in the limit of strong LPM 
suppression takes the form 



dl _ a s {A-Ax + 2x 2 ) I 2nC 3 {x) 



dxdz 3ir y E q x 3 (l — x) 

Ignoring the contributions to the energy loss from the two narrow regions near 
x ~ and x ~ 1, in which ( |4.81| ) is not valid, one finds the energy loss per unit 
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length 

« l.la 3 y/nC 3 (0)E q , (4.82) 

where to logarithmic accuracy p e jj ~ [a^n£^ 9 x(l — x)]~ 1//4 is taken. Note that 
as in QED, the elimination of the infrared divergence for strong suppression is 
a direct consequence of the medium modification of the gluon formation length 
which makes the typical transverse distances much smaller than l/m 9i9 . 

The medium modification of the formation length plays an important role in 
the case of gluon emission by a quark produced inside a medium. In this case the 
finite-size effects become important and suppress the radiation rate (cf. ( |3.47| )). 
This finite-size suppression leads to the 1? dependence of the quark energy loss 
for a high energy quark ( 3.49| ). One obtains 

AE q ~ a s C 3 (0)nL 2 . (4.83) 

This regime takes place as long as L < [Eq/nC^iO)) 1 ' 2 . Then it transforms into 
the AE q oc L behavior given by ( [4.82 ). More detailed discussions and numerical 



estimates are given in [31, 32]. The g — > gg transition can be evaluated in an 



analogous way. The ggg system can be in antisymmetric (F) and symmetric (D) 
color states. However, the two-gluon Pomeron exchanges do not generate the 
F <-> D transitions. This allows one to express the emission probability through 
the Green function for F state. 

4-4 Comparison with the BDMPS approach 

We conclude this section with a comment on the connection between the path 
integral approach with the approach discussed in section 3. 

Let us consider the case of a parton entering the medium from outside. The 
equivalence of the two approaches may be established using ( 4.72j )-( |4~74" ) together 
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with ( 4,7§| ), in the zero mass case as assumed in BDMPS pi]] . As it was men- 



tioned, the "frozen-size" expression (4.78) corresponds to the factorization con 



tribution, neglected in the BDMPS approach, on the ground of its weak medium 
dependence. Rewriting ( 4.72| ) as 



dl dl abs dP r ,dP r dI BH s 

+ -—-(- — ), (4.84) 



dx dx dx dx dx 
and ignoring the second term, one can show that 

dl dl abs 



(4.85) 



dx dx 

together with identifying in ( 4.74| ) the product ^ b ^aabc and 3> with the amplitudes 



/Bom and /, respectively, which are discussed in section 3. 

For the case of a parton produced inside a medium, say at z = 0, in ( 4.70 ) 



2; 



0, one should subtract from the right hand side of ( |4.72| ) the contribution 
corresponding to the configurations with z\ < and Z2 > in ( [4.70 ). The addi- 



tional term corresponds to the additional contribution in the BDMPS approach 
due to the hard scattering in the medium. In this case the "frozen-size" ex- 
pression is medium independent and it may be obtained by taking the limit of 
dl abs /dx when u> — > oo. Further details concerning this comparison can be found 



in 1 25] 



5 RADIATIVE ENERGY LOSS IN AN EXPANDING QCD 
PLASMA 

In the previous sections we have discussed the suppression of gluon radiation due 
to multiple scatterings of energetic partons propagating through dense matter 
with properties constant in time. 

Here we consider the case of a parton, of high energy E, traversing an expanding 
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hot QCD medium. We concentrate on the induced gluon radiation, the resulting 
energy loss of a quark, and its relation to jet broadening 37, 38]. 



Let us imagine the medium to be a quark-gluon plasma produced in a relativis- 
tic central AA collision, which occurs at (proper) time t = 0. We have in mind the 
realistic situation where the quark is produced by a hard scattering in the (not 
yet thermalized) medium, and at time to it enters the homogeneous plasma at 
high temperature To, which expands longitudinally with respect to the collision 
axis. Consider to to be the thermalization time, and for most of the results the 
limit to — > may be taken with impunity. The quark, for simplicity, is assumed 
to propagate in the transverse direction with vanishing longitudinal momentum, 
i.e. at rapidity y = 0, such that its energy is equal to its transverse momentum. 
On its way through the plasma the quark hits layers of matter which are cooled 
down due to the longitudinal expansion. It is assumed that the plasma lives long 
enough so that the quark is able to propagate on a given distance L within the 
quark-gluon phase of matter. 

As a consequence of the medium expansion the parton propagation in the 
transverse direction, z, is affected by the position-dependent density of the plasma 
p(z) and the parton cross section da/d 2 q±(q±, z). Therefore the screening mass 
H and the mean free path A depend on z. When the properties of the expanding 
plasma are described by the hydrodynamical model proposed by Bjorken [pgf] , 
one has the scaling law 

T 3 t a = const, (5.86) 

where the (proper) time t at rapidity y = coincides with the distance z on which 
the quark has propagated through the plasma. The power a, approximated in 
the following by a constant, may take values between and 1 for an ideal fluid. 
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Correspondingly, the transport coefficient q(t) defined as 

da M 2 (0~ 



«)M0/fttf|--* (5.87) 



<Pq± \{t) 

becomes time-dependent and satisfies 



m = ?(to) ( t ) > ( 5 - 88 ) 



due to ( gjHp . 

As a result [|37j the radiative energy loss AE for the quark (produced in the 
medium) traversing an expanding medium is 



AE = -^^q(L)L*. (5.89) 
2 — a 4 



In the high temperature phase of QCD matter |1C] 

l-a = 0{a 2 s {T)). (5.90) 

The coefficient q(L) = q(T(L)) has to be evaluated at the temperature T(L) 
the quark finally "feels" after having passed the distance L through the medium, 
which during this propagation cools down to T(L). One may, however, notice 
that the limit a = 1 for an expanding ideal relativistic plasma can be taken. 
In this limit the maximal loss is achieved. It is bigger by a factor 2 than the 
corresponding static case at fixed temperature T{L). 

So far we have discussed the result for the case for E > E cr (L), actually taking 



E — ► oo. In [38 1 the approach of the quark's energy loss AE to this limit is 
studied numerically as a function of the quark energy E. For instance, with 
L = 6 fm, one finds (almost) energy independence on E, when E > 100 GeV 



~ E cr , as given by fl5.89| ). 

In summary one expects indeed that the energy loss in an expanding medium 
be larger than in the static case taken at the final temperature, since the parton 
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passes through hotter layers during the early phase of the expansion. Perhaps 
the surprising feature is that there is no dependence of the enhancement factor 
on the initial temperature To. This result has to be associated to the coherence 
pattern of the medium induced radiation. Gluons contributing to the energy loss 
require finite time for their emission, and therefore effects of the early stages of 
the quark-gluon plasma expansion are reduced. 

6 INDUCED ENERGY LOSS OF A HARD QUARK JET IN A 
FINITE CONE 

Let us consider a typical calorimetric measurement of hard jets produced in heavy 



ion collisions [40|. The consequence of a large energy loss is the attenuation of the 
spectrum usually denoted as jet quenching. It is necessary to study the angular 
distribution of radiated gluons in order to give quantitative predictions for the 
energy lost by a jet traversing hot matter. Only the gluons which are radiated 
outside the cone defining the jet contribute to the energy loss. 



In |34|, 41, [42| the calculation of the angular distribution is discussed for a 
hard jet produced in the medium. Here we have in mind a hard quark jet of high 
energy E produced by a hard scattering in a dense QCD medium and propagating 



through it over a distance L. Following [41| we concentrate on the integrated loss 



outside an angular cone of opening angle 6 cone (Figure [|), 

- AE(e conc ) = l duj -r-r-r a de - ( 6 - 91 ) 

Jo Je couc duidzdv 
In the following we consider the normalized loss by defining the ratio 

R(e c0 n C ) = AJ ^° ne) - (6.92) 
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This ratio R(6 C one) turns out to depend on one single dimensionless variable 

R = R(c(L)6 cone ), (6.93) 

where 

c 2 (L) = ^-q(L/2f. (6.94) 

The "scaling behaviour" of R means that the medium and size dependence is 
universally contained in the function c(L), which is a function of the transport 
coefficient q and of the length L, as defined by ( |6.94| ). In Figure 10, we show the 



variation of R with 9 CO ne- The ratio R(6 C one) is also universal in the sense that it 
is the same for an energetic quark as well as for a gluon jet. The fact that 8 cone 



scales as l/c(L) may be understood from the following physical argument [22]: 
the radiative energy loss of a quark jet is dominated by gluons having uj ~ qL 2 . 
The angle that the emitted gluon makes with the quark is 6 ~ kj_/uj, and k"j_ ~ qL 
so that the typical gluon angle will be 9 2 ~ 1/qL 3 . 

So far we have discussed the medium-induced energy loss. Concerning the total 
energy loss of a jet of a given cone size it is important to take into account the 



medium independent part, which for a quark jet in a cone may be estimated [41], 



- AE^(9 conc ) ~ In (^) , (6.95) 

3 7T Vc'cone/ 

using a constant a s (6> max is taken 0(n/2)). 

7 PHENOMENOLOGICAL IMPLICATIONS 

The parameter controlling the magnitude of the energy loss is q. Estimates can 
be provided for its value, allowing us to give orders of magnitude for the radiative 
induced energy loss. The following numbers are estimates for a quark jet produced 
in matter. 
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For hot matter taking T = 250 MeV, 4- ~ 1 GeV/fm 2 taken from perturbative 



estimates at finite T, a typical value for v ~ 2.5, we find <f ~ 0.1 GeV 3 24]. With 
a s = |, this leads for the total induced energy loss to 



L ^ 2 



A£ « 60 GeV ( — ) (7.96) 



In [24 1 it is shown that for cold nuclear matter it is possible to relate q to the 
gluon structure function G evaluated at an average scale /J- 2 j^, actually 

27r 2 (T 

?~-^p[iG(4 (7.97) 

Taking the nuclear density p ~ 0.16 fm -3 , a s = |, xG ~ 1 for x < 0.1, it is found 
that 

L - 2 



- A£as4GeV lw^J ■ (7 - 98> 

These values do suggest that hot matter may be effective in stimulating significant 
radiative energy loss of high energy partons. As discussed in section 5 the energy 
loss is larger in an expanding hot medium than in the corresponding static one. 

Next we turn to the medium-induced — AE(9 cone ) for energetic jets. We may 
use the estimates above to give orders of magnitude for c(L) in the case of a 
hot/cold medium : 

c{L) hot ~ 40 (L/10 fm) 3 / 2 . 
A much smaller value is found in the cold nuclear matter case : 

c(L) cold ~ 10 (L/10 fm) 3 / 2 . 

As expected from the fact that 7?(^ cone ) depends universally on c(L)9 cone , Figure 



10 shows that the jets are more collimated in the hot medium than in the cold 
one. The loss is, however, still appreciably large even for cone sizes of order 
ft ~ 30° 
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Again, keeping in mind that the estimates are based on the leading logarithm 



approximation, we show in Figure 11 the variation of AE(9 cone ) with 9 cone of 
the medium-induced (for a hot medium with T = 250 MeV) and the medium 
independent energy losses. 

Let us now give a few representative examples of phenomena sensitive to par- 
ton energy loss in dense matter. Available experimental results are, it seems, 
essentially instructive for future measurements at higher energies. 

As a projectile traverses dense nuclear matter the width of the transverse mo- 
mentum distribution of partons may increase. pA and AA scattering allow to 
study parton pj_ -broadening of initial quarks in the Drell-Yan process of lep- 
ton pair production, and of gluons in J/^f production, respectively (see e.g. 
[ }43| , f44[). Recently, the analysis of J/\& data shows indeed p± -broadening of 
the intrinsic gluon distribution, which when translated into the ratio q/p results 



into q/p = 9.4 ± 0.7 [45|. Neglecting final state effects this should be compared 
with q/p ~ 7.4 [xG(x)\. 

The above quoted processes also contain information on the energy loss in the 



initial state due to matter effects [46, 47 j. In the Drell-Yan process the observed 
energy loss — AE of the incident quarks is indeed compatible with the estimate 
given in (|7.98| ), including the L 2 dependence, as measured and analysed in |48[| . 

Large p± particle and jet spectra and production rates in high-energy collisions 
are especially sensitive to a finite energy loss, when the partons are propagating 
through long-lived high density media before hadronization. Under extreme con- 
ditions the jets may even be "extincted" Q. However, hadron spectra from 
present experiments of pp, pA and AA collisions, mainly from CERN-SPS, do 
not show any strong evidence of suppression [E^, |50p. This observation which is 
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obscured by large theoretical uncertainties may indicate that at present energies 
the life time of the dense partonic matter may be shorter than the mean free 
path of the propagating partons. Significant jet quenching should become clearly 
observable in AA collisions at RHIC and higher energies, even for transverse mo- 
menta as low as p± > 3 GeV: the magnitude of the predicted jet quenching is 
commented upon in ^] . Suppression of hadronic p± distributions in the case 
of a thin plasma, therefore due to only a small number (< 3) of scatterings, is 



analysed in 53]. 



Jet quenching for very high energy jets is also discussed in [40, 55]. In 
particular the ratio of monojet versus dijets observed in ultra-relativistic heavy 
ion collisions is predicted. 

A further interesting proposal to study the modification of jet fragmentation 



due to energy loss is proposed in [56, 57]. Noting that photons are essentially 
not affected by hadronic media, the conjecture is to measure the charged particle 
p± distribution in the opposite transverse direction of a tagged photon, i.e. in 
7 + jet events of high energy heavy ion collisions. With increased luminosity this 
may be even possible at RHIC energies. 

Valuable and important information about dense hadronic matter produced 
in collisions is provided by dileptons, either from Drell-Yan processes or from 



final heavy meson decays [58]. In this context Shuryak j59[ pointed out the 
importance of the energy loss of charm (bottom) quarks due to their interactions 
in the medium. In the extreme case they may be even stopped in dense matter. 
Assuming e.g. that the charmed mesons D and D* take all the charm quark 
momentum in the fragmentation process the final leptons from the semileptonic 
decay populate the invariant dilepton mass spectrum at masses below 1 — 2 GeV 
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(4 — 5 GeV from bottom decays): as a result dilepton spectra in AA collisions 
for invariant masses above 2 GeV are not dominated by correlated semileptonic 
charm and bottom decays. This expected strong suppression due to energy loss 
is confirmed in further detailed (Monte Carlo) studies in [60, [H], 62]. 



8 OUTLOOK 

In this review we have described the more recent results related to energy loss 
and pj_-broadening of a high energy quark or gluon (jet) traversing QCD media. 
Phenomenological implications for measurements in cold as well in hot (QGP) 
matter have been discussed. The orders of magnitude found for the energy lost by 
an energetic jet in hot deconfined matter indicate the interest of the corresponding 
measurements as specific signals. 

A couple of important open questions triggered by the coherent character of 
the induced energy loss remain open. One is related to the formulation of a 
transport model (Monte Carlo) which correctly simulates the interference pattern 



of gluon radiation induced by multiple scattering [18, 22]. It is indeed crucial when 



calculating rates for processes leading to thermal and chemical equilibration of 



partons to include medium effects [63] - J6S|] . In the same context let us mention 
the influence of the LPM effect on the production of dileptons and real photons 
produced in a QGP, or in a hadron gas |ffo| . In any case the partons are not very 
energetic, since their energies are of the order of the plasma temperature. This 
forbids to use the asymptotic treatment discussed in the above. 

We already mentioned that the above discussed numerical estimates give only 
orders of magnitude in particular since they are obtained in leading order in the 
QCD coupling. The main aim of the present investigations is therefore to encour- 



39 

age experimentalists at RHIC, and later at LHC, to carefully explore heavy quark 
production and especially jet phenomena in ultra-relativistic heavy ion collisions 
||71[|. High p± nuclear physics may become an exciting new frontier at these col- 



liders [52], because of the possible jet "extinction" or crucial modifications of the 
spectra 0. The best guiding example comes from the long-term study of hard 
jets in hadron-hadron scattering, starting from the first evidence at the CERN- 
ISR until the analysis of jet cross sections up to transverse momenta of p± ~ 500 
GeV at CDF and D0, which has been successfully carried out within an active 



interplay between experiments and perturbative QCD [72]. 

Medium effects will continue to attract increasing attention, and we hope that 
finally the described suppression mechanism will be demonstrated by future ex- 
periments in a convincing manner. More detailed treatments and improvements 
are certainly necessary in view of completing this program. It constitutes an im- 



portant chapter of what has been recently designated [73] as the "health report" 
of QCD. 
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Figure 1: Self-energy diagrams contributing to the collisional energy loss: (a) in 
HTL-resummed perturbation theory for soft exchanged momentum and (b) in 
fixed leading order for hard exchanged momentum. 
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Figure 2: Collisional energy loss of a charm quark as a function of its momen- 
tum. The quark propagates through an out-of-chemical equilibrium plasma with 
fugacities X g = l,X q = (solid curve) fl~6[| . The dashed curve is the equilibrium 
result of []n|, the dotted curve shows the original prediction by Bjorken . 
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MEQIUM 



Figure 3: Contribution to the induced gluon spectrum by interference between 
the amplitude / and the Born amplitude fsom ■ 
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Figure 4: Energy loss as a function of the medium size L. 
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Figure 6: The diagrams representing the Born terms for the amplitude. 
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Figure 7: Graphs describing the gluon emission. 
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Figure 9: Example of a hard process producing a quark jet. The gluon is emitted 
outside the cone with angle CO n e . 
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Figure 10: Medium induced (normalized) energy loss distribution as a function 
of cone angle con e f° r hot (T = 250 MeV) (solid curve) and cold matter (dashed 
curve) at fixed length L = 10 fm. 
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Figure 11: Energy loss in a hot medium, T = 250 MeV, as a function of 9 cone . 
The dashed curve represents the medium independent piece for E = 250 GeV. 
L = 10 fm. 



